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A self-consistent theory involving Maxwell equations and a density-matrix 
linear-response theory is solved for an electromagnetically-coupled doped 
graphene micro-ribbon array and a quantum-well electron gas sitting at an 
interface between a half-space of air and another half-space of a doped semi- 
conductor substrate which supports a surface-plasmon mode in our system. 
The coupling between a spatially-modulated total electromagnetic field and 
the electron dynamics in a Dirac-cone of a graphene ribbon, as well as the 
coupling of the far-field specular and near-field higher-order diffraction modes, 
are included in the derived electron optical-response function. Full analytical 
expressions are obtained with non-locality for the optical-response functions of 
a two-dimensional electron gas and a graphene layer with an induced bandgap, 
and are employed in our numerical calculations beyond the long-wavelength 
limit (Drude model). Both the near-field transmissivity and reflectivity spec- 
tra, as well as their dependence on different configurations of our system and on 
the array period, ribbon width, graphene chemical potential of quantum-well 
electron gas and bandgap in graphene, are studied. Moreover, the transmit- 
ted E-field intensity distribution is calculated to demonstrate its connection 
to the mixing of specular and diffraction modes of the total electromagnetic 
field. An externally-tunable electromagnetic coupling among the surface, con- 
ventional electron-gas and massless graphene intraband plasmon excitations 
is discovered and explained. Furthermore, a comparison is made between the 
dependence of the graphene-plasmon energy on the ribbon width and chemical 
potential in this paper and the recent experimental observation given by Ju, et 
ai, [Nature Nanotechnology, 6, 630 (2011)] for a graphene micro-ribbon array 
in the terahertz-frequency range. 
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1. Introduction 

There has been a high-level of research and device interests on both the electronic and 
optical properties of two-dimensional (2D) graphene material [IH5] since the report on the 
first successful isolation of single graphene layers as well as the related transport and Ra- 
man experiments performed for such layers [6] . The most distinctive difference between a 
graphene sheet and a conventional 2D electron gas (EG) in a semiconductor quantum well 
(QW) is the electronic band structure, where the energy dispersions of electrons and holes 
in the former are linear in momentum space but they become quadratic for the latter. As 
a result, electrons and holes in graphene behave like massless Dirac fermions and show a 
lot of unexpected physical phenomena in electron transport and optical response, includ- 
ing anomalous quantum Hall effect [7J IB], bare and dressed state Klein tunneling [9l-fT2"] and 
plasmon excitation pT3HT5| . a universal absorption constant [16j U7] . tunable intraband [18] 
and interband [THIED] optical transitions, broadband p-polarization effect [21], photoexcited 
hot-carrier thermalization (22] and transport [23], electrically and magnetically tunable band 
structure for ballistic transport [21] , field-enhanced mobility in a doped graphene nanorib- 
bon [22] and electron-energy loss in gapped graphene layers [2S] ■ 

It has been found that most of the unusual electronic properties in graphene can be 
very well explained by single-particle excitation of electrons. For diffusion-limited electron 
transport in doped graphene, the Kubo linear- response theory [2Zj, as well as the Hartree- 
Fock theory in combination with the self-consistent Born approximation [28j , were used. 
Moreover, the semiclassical Boltzman theory was applied in both the linear [2H] and non- 
linear [25] response regimes. For the collective excitation of electrons in graphene, previous 
studies found it has an important role in the dynamical screening of the electron-electron 
interaction [1314151 l26| 130] . However, the EM response of graphene materials has received 
relatively less attention, especially for low-energy intraband optical transions [T8l l3T| 132] . 

For any charged-particle systems, there always exists an electronic coupling based on 
the intrinsic electron-electron interaction in addition to the system response to an external 
electromagnetic field or an external electron (ion) beam. This electronic coupling can be 
either a classical direct coulomb interaction between any two charged particles with no 
requirement for an overlap of their wave functions, or quantum-mechanical tunneling and 
exchange coulomb interaction with a tail overlap of their wave functions [33] . For the same 
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charged-particle systems, there exists another electromagnetic coupling based on the light- 
electron interaction [31] . This electromagnetic coupling results from an induced optical- 
polarization field in the system due to dipole coupling to an incident light wave. The total 
electromagnetic field, including the induced optical-polarization field, should satisfy the 
macroscopic Maxwell equations [35] . However, the induced optical-polarization field in the 
system is determined by the microscopic density-matrix equations [36J. Therefore, a self- 
consistent theory is required for studying the electromagnetic response of such a charged 
system to incident light. 

Physically, a statistical ensemble average of dipole moments in an electronic system with 
respect to all quantum states will give rise to an induced optical-polarization field [37]. This 
provides a connection between the microscopic density-matrix equations and the macroscopic 
Maxwell equations |36j . The former describes the full electron dynamics excited by a total 
electromagnetic field, including optical pumping, energy relaxation, optical coherence and 
its dephasing, and carrier scattering with impurities, lattice and other charged carriers p5] 
HO]. When the incident field is weak, a self-consistent linear-response theory to the total 
electromagnetic field (not just the incident light) can be applied to an electronic system as 
a leading-order approximation. 

Experimentally, the inelastic light scattering [41J and electron-energy-loss spec- 
troscopy [12] were widely used for probing the bulk and surface electronic coupling, re- 
spectively, in a charged-particle system. On the other hand, the elastic light scattering [13] 
and the optical absorption spectroscopy [H] were extensively applied for studying the sur- 
face and thin-film electromagnetic coupling in either a conductive or a dielectric system. 
For a smooth conductive surface, only electronic excitations in the long-wavelength limit 
can be coupled to an incident light wave, where a Drude-type model is usually adopted 
in spectroscopy analysis. For a patterned conductive surface in the deep sub-wavelength 
regime, on the other hand, surface light scattering can be used to explore surface electronic 
excitations in the system beyond the long- wavelength limit, where a wavenumber-dependent 
(nonlocal) optical response of electrons, which is crucial for the near-field distribution and 
the Landau damping of plasmons [33], must be employed to analyze the light reflectivity and 
transmissivity spectra. 

In this paper, we consider a model in which a two-dimensional electron gas (2DEG) and 
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a doped linear graphene micro-ribbon array (electronically isolated from each other) are put 
on the surface of a doped semi-infinite semiconductor substrate. Therefore, we expect to see 
the excitations of the localized surface plasmon and the conventional 2DEG plasmon, as well 
as the 2D Dirac-fermion plasmon excitation, in such a system in the presence of an incident 
electromagnetic wave. More interestingly, we expect there to exist strong electromagnetic 
couplings among these three collective excitations in the deep sub-wavelength regime, i.e. 
with an array period much less than the incident wavelength, and these couplings can be 
externally tuned by the doping densities in a graphene layer or in a QW, as well as by the 
driven 2DEG in the QW or by the induced bandgap in the graphene. The existence of 
an bandgap in the graphene not only modifies the Landau damping of plasmons |15j. but 
also introduces the interplay between Dirac-like and Schrodinger-like electrons in the elec- 
tromagnetic response. To demonstrate the predicted strong and tunable electromagnetic 
coupling among three plasmon excitations in our model system, numerical calculations and 
their comparisons are presented for both the transmissivity and reflectivity spectra with dif- 
ferent configurations of our model system and various array periods, ribbon widths, chemical 
potentials of graphene, and induced bandgaps in graphene. 

The rest of the paper is organized as follows. In Sec.[2j we derive the self-consistent 
Maxwell equations and the density-matrix linear-response theory within the random-phase 
approximation for our model system which consists of a combined 2DEG in an InAs QW 
and a linear graphene micro-ribbon array doped by electrons on the surface of a semi-infinite 
n-doped GaAs substrate. This self-consistent theory, along with the detailed expressions for 
a linear matrix equation in Appendix A, is then applied to study the transmissivity and 
reflectivity spectra of an incident plane wave for various structural-geometry and material 
parameters. In Sec.[3j numerical results, including the p-polarized far- and near-field trans- 
missivity spectra as well as the transmitted p-polarized E-field intensity distribution, are 
presented to demonstrate and explain the existence of an externally-tunable electromag- 
netic coupling among the localized surface plasmon, 2DEG plasmon and massless graphene 
plasmon excitations in our system. The conclusions drawn from these results are briefly 
summarized in Sec.|U 
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2. Model and Theory 

It is shown in Fig. [I] that the half-space z < is air with the refractive index n a — 1 
whereas the half-space z > is occupied by a doped semi-infinite GaAs substrate with 
a dynamical dielectric function e s (q, oS) which supports a localized surface-plasmon (SP) 
mode. At the interface z — 0, there are a periodic electron-doped graphene micro-ribbon 
array (GMRA) and an InAs quantum well (QW) for a 2DEG. Optically, the GMRA and 
QW can be considered to be in one plane. The reason for this is that the sheet thickness 
is much smaller than the wavelength and the decay length of the impinging light. They 
can, however, still be electrically separated from each other by an energy barrier [3U . 
A plane-wave electromagnetic (EM) field incident from z < is diffracted by the GMRA 
into both the reflection (z < 0) and transmission (z > 0) regions. The diffraction modes 
of the EM field are modulated by the induced optical polarization at the interface with the 
GMRA, which will be electromagnetically coupled to SP and 2DEG plasmons. 

We assume that the GMRA is periodic in the x direction, and seek the solutions based 
on a plane-wave expansion for a single frequency. Setting q y = for simplicity and denoting 
the transverse electric component of the EM field below the interface by E< and above by 
E>, we obtain [Mffl 
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Here, n = 0, ±1, ±2, • • • label all reciprocal lattice vectors nG = n(27r/d), d is the period 
of the GMRA, q n = q x + nG, q x = (co/c) n a sin^, n = (u/c) n a cos 9i, n a is the refractive 
index of air, 9{ is the angle of incidence, and u the angular frequency of the incident EM 



field. In addition, (3i ;Tl and f}% >n i n Eqs. ([I]) and (|2j are determined by the following dispersion 
relations for transverse EM fields 
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where c is the speed of light in vacuum, e s (q x , (3, cu) s is the dielectric function of the doped 
semiconductor substrate, and we require that Im(/3i. n ), Im(/?2,n) > so as to ensure a finite 
value for the reflected or transmitted EM field at z = ±oo. 

Formally, if we keep only the term with n = in the summations in Eqs. ([TJ and ([2]), 
the GMRA in our model will simply change to a graphene sheet. For this situation, after 
removing the doping in the GaAs substrate and the 2DEG in the InAs QW, the formalism 
in this paper will reduce to the previous ones [3TJ E2] for a transverse EM response from a 
single graphene sheet. 

In the long-wavelength limit, the transverse dielectric function e s (q x ,(3, u) introduced in 
Eq. ^ for the doped semiconductor substrate is found to be 



e s (g x ,/3; u) 



1 - 



n 2 



(4) 



u(u + 270) 

where Q p \ = ^/n^e 2 / e ebm* is the bulk plasma frequency, n 3D is the doping concentration 
in GaAs, /ryo is the homogeneous level broadening, m* is the effective mass of electrons in 
the substrate and e& is the relative dielectric constant of the semiconductor substrate. 

Once the electric field E is obtained by using the Maxwell equations, the magnetic field 
H for non-magnetic materials can be obtained from 
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where jiQ is the vacuum permeability. For q y = and z < 0, Equation (|5| explicitly leads to 
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H<(x, z\ oo) 



UfJLQ 
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In a similar way, the magnetic field H> in the region of z > can also be obtained from 
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Electronically, the GMRA and 2DEG can be separated by an undoped wide-gap semi- 
conductor barrier layer. Optically, however, the GMRA and 2DEG can both be regarded 
as two-dimensional sheets since their thickness are much less than the EM-field wavelength 
considered. Therefore, the four independent boundary conditions for transverse EM fields 
at the interface z = can be written as I 
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where P s = (P x , Pf) is the induced sheet optical-polarization field (related to the po- 
larization current) while ot s = (a x , a y s ) is the induced sheet conduction-current density. 
Furthermore, we can write P s (x| oo) as [35j [36] 
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w)/e 




w)/e 



XiiQn, u) 



C x (q n \ co) 
C y (q n \ oo) 



S 



+ C Xa(gn', w)sinc[(n' 



Cy(?n'| W) 



(17) 



In Eq. (17), £ = W/d, W is the width of a graphene micro-ribbon, sinc(x) = sinx/x, and 
Xi(ga;, cu) and XziQx, w ) are the 2DEG and the graphene-ribbon polarizabilities, respectively. 
The multi-ribbon effect represented by the summation over n', as well as the mode- mixing 



effect with n ^ n', are both taken into account by the second term of Eq. (17). Moreover, 



the coupling between a spatially-dependent total EM field and the electron dynamics in an 
energy band, due to the non-locality in an optical response of electrons, is also included in 



Eq.(17) 



Similarly, using Ohm's law we can write the sheet conduction-current density as 



a a (x\u)= oc s (q n \ cj) e iqnX , 
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and ai(q x , u) and ^(fo, u) are the real 2DEG and graphene-ribbon optical conductivities, 
respectively. 

For an incoming p-polarized EM field, we get A x (q x \ uj) = T]H / (weon„) and A y (q x \ u) = 
0, where H is the magnetic component and €q is the permittivity of free space. 



Combining the sheet optical polarization in Eq. (16) and the sheet-current density in 



Eq. (18), we can write the optical-response function for j = 1, 2 in a compact form given by 



Xj(q, w) = Xj(q, + 



uie 



(20) 



where q = yjq\ + q\. 

We notice that there exists an anisotropy in the polarizability of a graphene nanoribbon 
along the longitudinal and transverse directions (13] , respectively, and it can be attributed to 
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a subband quantization due to finite-size effect. However, this anisotropy can be neglected 
for the width of ribbons in the micrometer range. In other words, we adopt in this paper the 
isotropy of the micro-ribbon array instead of anisotropy of the nanoribbon array. Within the 
random-phase approximation, for graphene ribbons in the micrometer range we calculate 
the nonlocal (q ^ 0) optical- response function Xziq, w ) as [E9 EH] 

e 2 

x 2 (g, u) = — 2 n 2 (g, u) , (21) 

e o<? 

where the non-retarded limit with uj ^> qc has been used and IT 2 (g, u) is defined as 



nl \ ^ \ \i ii -iar I l i \l 2 /o(^m,k) /o(^n2,k+q) / nr) \ 
2 {q, u = — > <?ii,k e q n 2 ,k + q — f— - , 22 

A ^ 1 k+q - e ni k - + «0+) 

ni,ri2, k ' M ' 

^4 is the surface area of a graphene micro-ribbon, fo(x) is the Fermi-Dirac function for 
thermal-equilibrium electrons, e ni k is the kinetic energy of electrons, and n\, n 2 = ±1 rep- 
resent the valence (— ) and conduction (+) bands of Dirac cones in a graphene layer. 

The Dirac-cone model only holds for crystal momentum that is not too large. In fact, the 
energy band structure ceases to be a cone above about 2.5 eV. For our medium-high-density 
sample used in this paper, the Fermi momentum is much less than the momentum separation 
between valleys at K and K' . The Coulomb interaction between electrons in different valleys 
will be significantly suppressed due to this large momentum separation. Additionally, the 
electron-phonon mediated inter-valley scattering also requires a temperature higher than 
room temperature. Since we assume very low temperature for electrons in our model, this 
inter- valley scattering may be justifiably neglected. 

The 2D dielectric function n(q, u) for graphene micro- ribbons can be obtained by [16] 
n(q, oj) = Kb + qX2{q, w)/2 in the non-retarded limit with ^ being the dielectric constant of 
graphene materials. In addition, the optical conductivity cr 2 (g, u) of graphene ribbons can be 



calculated from Eq. (20) within the linear-response theory, i.e. cx 2 (g, u>) = we Im[x 2 (g, lu)]. 



In a similar way, <Ji(g, oo) can be calculated from Eq. (41) below. Note that qX2(q, along 



with Eq. (21 ) correspond to 2D rather than 3D Fourier transform of the Coulomb interaction. 



After a lengthy calculation [13] . from Eq. (22) we get an analytical expression for a gapped 
graphene layer at T = K as follows: 
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n / \ _2^2 q 2 

x {i [G>(xi _) - G > (x li+ )] Q 1< (a; 2> _) + [G<(xi_) + iG > (x 1)+ )] Q 2< (x 2 _, x 2 ,+) 

+ [G<(.x li+ ) + G<(xi,_)] Q 3< K-) + [G<(a;i,_) - G<(x lj+ )] Q 4< (a; 2 , + ) 

+ [G?>(ar li+ ) -G>(xi_)] Qi > (x 2 ,- > a?s) + [G>(aJi,+) +zG<(xi_)] Q 2> (x 2 _, x 2 ,+) 

+ [G>(x lj+ ) - G>(-x L _) - ztt[2 - arjj]] Q 3> (a;2,+) 

+ [C>(-xi_) + G>(xi i+ ) - i7r[2 - Xg]] Q 4> (a; 2 -, x 3 ) 

+ [GoK+)-G'oK_)]Q 5> (x3)} , (23) 

where /i 2 = a/ (hvFk 2 F) 2 + {Eq/2) 2 — Eg/2 is the chemical potential of electrons in a 
graphene layer with respect to zero energy at k — 0, t>i? is the Fermi velocity of graphene, 
the kinetic energy of electrons in valence and conduction bands are 



e ± , w = ±JVv 2 F k 2 +E G /4, (24) 



and Eg is the induced bandgap. From Eq. (24) we know that a finite effective mass Eg/2v 



of electrons will be created for a gapped graphene layer close to the edge of a gap with 
k < E G /hv F . 

The coulomb interaction between electrons will induce a self-energy, leading to an energy- 
dependent Fermi velocity beyond the random-phase approximation j47j. However, this 
energy-dependent behavior becomes significant only around the Dirac point (48] . For our 
samples in this paper with a very high doping concentration, we expect this effect is min- 
imized. Therefore, we have neglected the coulomb renormalization to the graphene Fermi 



velocity in Eq. (24). Correction to the self-energy comes from correlation effects. However, 
this correction becomes significant only at low densities, as it occur, for example, when there 
is Wigner crystallization. However, the sample considered in this paper is of medium-high- 
density concentration, where correlation effects may be neglected. Since we recognize that 
the limitations of the model is a poor argument when it comes to the experimental data, we 
included in Eq. (24) an effective mass ~ Eg term into Dirac fermion picture. We note that 
the induced bandgap may be attributed to either photon dressed-state effects by circularly 
polarized light or to a substrate. There exist two main sources of disorder in graphene rib- 
bons, i.e., edge roughness and charged impurities. Since wide micro- ribbons are considered 
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in this paper, we can neglect the edge roughness. This leaves the charged impurities as the 
dominant scattering mechanism. The existence of charged impurities is also able to induce a 
self-energy, leading to a deformed density of states around the Dirac point of graphene [49J . 
However, for our medium-high-doping samples, the screening length of impurity potential 
is small, leading to a negligent effect on the density of states in graphene. Moreover, with 



the added gap in Eq. (24), the charged impurities create localized states only in the range 
between 0.92 E G and ^gISO]- Since our transmission is not substantially affected up to 
E G = 0.25 eV (see Fig. 11 in Sec. Ill) and the induced states by localized impurities are 
well above the terahertz regime of interest in this paper, we can neglect this effect here. 
For impurities, the interference or localization effect on the conductivity becomes important 
only for very high impurity concentrations. Since the sample discussed in this paper is only 
a medium-high-density one, localization is not expected to be a dominant effect in compar- 
ison with other effects studied in this work. The detailed discussion related to the screened 
disordered effect on the electron polarizability of graphene beyond the random-phase ap- 
proximation is, however, outside the scope of the current paper. 



Three functions introduced in Eq. (23) are defined as 



G<(x) = xyjxl - x 2 - (2 - x 2 ) cos" 1 (jpj , (25) 
G>(z) = x^jx 2 - xl - (2 - xl) cosh" 1 (j^J , (26) 



G (x) = xJx 2 -x 2 -(2- x 2 ) sinh- 1 —= . (27) 




Moreover, nine region functions used in Eq. (23) are defined by 
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Qi< O2,-) = <9(// 2 - x 2 ,_ - fiw) , (28) 
Q 2< (x 2 _, a?2,+) = 9(—fwj — fM 2 + x 2 ~) 0(hu + /x 2 — £2,-) 0(/x 2 + £2,+ - ^) , (29) 

Q 3< (s 2 ,-) = 0(-^2 + z 2j - - M , (30) 

Q4< (a?2,+) = + A*2 - ^2,+) 0(hv F q - frw) , (31) 

Qi > (x 2 _, £ 3 ) = 9(2k 2F - q) 6(hui - x 3 ) 9{\i 2 + x 2 _ - foj) , (32) 

Q2> (a?2,-, ^2,+) = 0(hu x 2 ,-) ^(a*2 + a; 2 ,+ - ^) , (33) 

Q3> (^2,+) = - A*2 - a?2,+) , (34) 

Q4> (#2 ,- , x 3 ) = 6*(g — 2£; 2F ) #(ftu; - sc 3 ) 0(/i 2 + ar 2 _ - TvJ) , (35) 

25> (x 3 ) = 9{hw - hv F q) 9(x 3 - frhj) , (36) 



where k 2F = {jx 2 + E G /2) 2 — (E G /2) 2 /hv F is the Fermi wave number and 9{x) is the unit 
step function. Finally, we have defined six variables Xq, £i,±, x 2) ± and x 3 through 



x = 

Xl,± 

x 2 ,± 
X3 = 



'1 + 



E 2 

G 



h 2 v 2 F q 2 — h 2 



2 ' 



2yU 2 ± ftjJJ 

Hv F q 



h 2 v 2 F (q±k 2F ) 2 + E 2 G /4 



h 2 v 2 F q 2 + E 2 G . 



(37) 

(38) 

(39) 
(40) 



Following a similar approach, within the random-phase approximation we get the nonlocal 
optical response function for the 2DEG inside an InAs Q W [IS] 



2p s e 2 m* s 
e h 2 k 1F q 3 



2z-C^ v / (z- U y 



C + ^{z + u) 2 



+ i 



£L V 1 ~ ( z ~ «) 2 " D + V 1 ~ ( z + u ? ] } 



(41) 



where p s is the 2DEG density, k\ F = ■^/2np s and m* are the Fermi wave number and the 

; in Eq. (41 ) for a driven 
(z+u)/\z+u\ and D + = (C+ = and £>+ = 1) 



effective mass of 2DEG. Additionally, we have defined the notations in Eq. (41) for a driven 
2DEG: u = m*u/hqk 1F , z = q/2k 1F , C + 
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for \z + u\ > 1 (\z + u\ < 1), and C- = (z - u)/\z - u\ and _D_ = (C_ = and _D_ = 1) 
for \z — u\ > 1 — u\ < 1). 



For the specular mode with n = 0, from Eqs. ( 12 )— (15 ) the boundary conditions for the 
EM fields at the interface yield 



B x (q x \u) - C x {q x \ u) = - A x (q x \u) 



(42) 



B y (q x \ u) - C y {q x \ u) = - A y (q x \ u) 



(43) 
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i\tU) 
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oo 
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n — — oo 



Xi(g«, w) + Cx2 (<?*•, w) — 4t(&| v) 



(45) 



The terms with n' ^ in Eqs. (44) and (45) reflect the coupling between the specular mode 



and higher-order diffraction modes. These generated higher-order diffracted EM fields can 
be treated as source terms in addition to the incident EM field in a perturbation picture. 



If we exclude all the terms with n' ^ in Eqs. (44) and (45), the solution of Eqs. (42)-(45) 



gives rise to a transverse EM response for a single graphene sheet by taking Xi(Qx, w ) = 0, 
C = 1 and e s (q x , p\o, w) = e 6 [HU E2]. 



For the diffraction modes with n = ±1, ±2, • • • , on the other hand, from Eqs. ( 12 )— ( 15 ) 
the boundary conditions for the EM fields lead to 



B x (q n \u) - C x {q n \uj) = 



(46) 



B y (q n \u) - C y (q n \uj) = 



(47) 
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ic 2 ^ 



— BJqJ u) - 



2.11 



uj- 



C y (q n \ w) 



C sinc[(n - n)ir(} xs(q n ', w) w) - Xi(?n, w) B tf (g„| w) = , (48) 



n =— oo 



n 2 a B x (q n \ w) - -5 — e s (q n , /3 2 , n , C^nl cu) 



oo 

-C sinc[(n' - n)7rC] x 2 (?n', w) B x (q n >\ u) - Xi(g n , u) B x (q n \ u) = , (49) 

n'=— oo 

where xi — an d Xi 7^ correspond to a GMRA only or a 2DEG plus a GMRA, respectively. 



The terms with n' = in Eqs. (48) and (49) reflect the coupling between different diffraction 



modes with the specular mode. Here, the terms with n' = 0, excited directly by the incident 
EM field, can be regarded as source terms for the generated higher-order diffracted EM fields 
with n = ±1, ±2, • • • in a perturbation picture. 



Once the solution of the matrix equation derived from Eqs. (42)-(49) (see Appendix A) 



is obtained, for p polarization we can express the intensities of EM field at the interface as 



\E x (x,z = 0| u)\ 2 = [A x (q x \ u)} 2 + ^ \B x (q n \uj)\ 2 + 2A x (q x \u)Re[B x (q x \uj)} 



+ B x (q n \co)[B x (q n ,\co)Ye^-^ Gx 

n, n' 

+ 2A x (q x \ u) Re ( ^ B x (q n \ u) e mGx j 



(50) 



\E y (x,z = 0\co)\ 2 = [A y (q x \u)} 2 + \B y (q n \co)\ 2 + 2A y (q x \co)Re[B y (q x \co)] 

n=— oo 

(n^n') 



+ J2 B y (q n \oj)[B y (q nl \uj)]*e^-^ Gx 

n, n' 

+ 2A y (q x \u)Re ( ^ B y (q n \ u) e inGx J 



\E z (x, z = 0| to 



2 00 2 

\ 2 = q f 2 [A x {q x \uj)} 2 + ]T t£^\BMu)\< 



(51) 
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o 2 ( n ^«') 



3 Ae^l W) Re I W)] + - 77T U B x (q n \ Co) [B x (q n >\ Co))* 

V i Pl.nVPl.n' ) 

n,n' ' ' 



X 



(n^O) 

S n-n>) Gx _ 2Ax (q x \co) Re | J2 ^ B x {q n \ co) e mGx ] , (52) 



n 



1,11 



where the interferences taking place in the x direction at the interface between the reflected 
EM field (n ^ n') and that between the incident and reflected EM fields (n ^ 0) are all 
included. 



Using Eqs. (50)-(52), we get the x-dependent square of the ratio lZ(x\ co) of the reflected 



to the incident E-field amplitude at the interface, given by 

| v z = 0| co)\ 2 + \E y {x, z = 0\co)\ 2 + \E' z (x, z = 0| co)\ 2 

TZ < [x\ co) = i , 53 

1 ' ' (l + Q 2 Jv 2 )Al(q x \u) + Al(q x \u;) 1 ' 

where lZ^(x \ co) contains both the near (<) and far (>) field contributions, and 

\E r x (x,z = 0\u)\ 2 = J2 B Mco) [B x {q n ,\uj)]* e^-^ Gx , (54) 

n, n' 

\E r y (x,z = 0\uj)\ 2 = J2By(qn\uj) [B y (q n ,\u)]*e^ n -^ Gx , (55) 



\E:(x,z = 0\co)\ 2 = £ a v B x {q n] co) [B x {q n ,\uW e^ G * . (56) 

— , Pl,n{Pl,n') 
n,n' 

Following the same approach, we get the x-dependent square of the ratio J-(x\co) of the 
transmitted to the incident E-field amplitude at the interface, given by 

= \Ej(x, z = 0| uo)\ 2 + \El(x, z = 0| co)\* + \Ej(x, z = 0| uo)\ 2 
[X{UJ) ~ (l + q 2 Jv 2 )AUq x \co) + Al(q x \co) ' 1 J 

where J r ^-{x \ co) also includes both the near (<) and far (>) field contributions, and 



\E t x (x,z = 0\co)\ 2 = J2cI(q n \oo) [C x {q n ,\ co)}* e^"'^ , (58) 



\El(x, z = 0| co)\ 2 = C y (q n \ u) [C y (q n ,\u)]* e^ n ~^ Gx , (59) 

n, n' 
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\El(x,z = 0| u)\ 2 = V - ^ a(gnk) [a(^|a;)]*e^-"') fe . (60) 

— , P2,n[P2,ri) 



n, tv 



After integrating Eqs. (53) and (57) with respect to x, we obtain the far-field reflectivity and 



transmissivity spectra of the sample as 



i f 

1Z(u) = lim — — / dx TZ > {x\ u) = } R n (q n \u) 9(n a uj - \q n \c) , (61) 



If 

F{uj)= lim — / dxJF > {x\uj) = V] F n (q n \Lu)9(^/7 s LU-\q n \c)6(e' s ) , (62) 

L x ->oo IL X J *— ' 



where e' s = Re(e s ), the interference terms between different diffraction modes are canceled 
out and 



R n (q n \ co) = lir ^ |W -" Nl ^" 1 ^ 1 Z^v^pi , (63 ) 



(l + ql/\Pl,n 


2 ) \B x (q n 


oo)\ 2 + \B y (q n 


u)\ 2 


(1 + ql/V 2 ) Al(q x \ 
(l + \q 2 /\^ n \ 2 ) \C x (q n 


w) + A 2 y (q x \to] 
\u)\ 2 + \C y (q n 


u)\ 2 



3. Numerical Results and Discussions 

In our numerical calculations, we assume p polarization for the incident EM-field and set 
the parameters in numerical calculations as follows: p s = 8 x 10 n cm -2 , m* = 0.024 mo 
with m being the free electron mass, E G = 0, fi 2 = 0.45 eV, vp = 10 6 m/sec, e& = 12, 
Mlpi = 10.14 meV, h^o = 5meV, ( = 0.5, d = 4 /im (deep sub-wavelength regime), 9i = 30° 
and q x = {to/c) n a sin 9\ with n a = 1. In this paper, we choose the color scale from blue 
(darker, minimum) to red (lighter, maximum). The change of these parameters will be 
directly indicated in the figure captions. In this paper, we only consider the low-energy 
intraband plasmon excitation in the terahertz-frequency range. 

To uncover the physical mechanism behind the reported graphene-plasmon peak shift with 
GMRA period [181 . we show a comparison in Fig.[2|for two calculated far-field transmissivity 



spectra J>(w), defined by Eq. (62), as a function of photon energy hw with d = 2/im (red 
solid curve) and d = 4 p,m (blue dashed curve). Here, only the graphene micro-ribbon array 
is included in our system. As d increases, the graphene plasmon peak at hu = 19.8 meV 
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for d = 2 /an shifts down to hw = 14.0 meV for d = 4/tm and this peak shift accurately 
satisfies the ~ 1/ \[d scaling relation experimentally observed [TS] . A series of kinks below 
the plasmon peak in Fig. [2] corresponds to the single-particle excitations at hw = hvp\qn\ (for 
\Qn\/k2F < 1) with q n « n(2ir/d) for different values of integer n. On the other hand, the 
shoulder appearing above the dominant peak comes from another severely-Landau-damped 
graphene-plasmon peak, which becomes more visible for d = 2 /im since the peak separation 
is increased. 

To see the many-body effect on the hybridized GMRA-QW plasmon excitation, Fig. [3] 
presents a comparison for J r p {w) as a function of hw with two values for the graphene chemical 
potential, i.e., /a 2 = 0.45 eV (red solid curve) and /x 2 = 0.9 eV (blue dashed curve). Here, the 
chemical potential /i 2 of graphene ribbons can be tuned by a gate voltage [18]. As expected, 
one of the two graphene-like plasmon peaks at hw = 19.5 meV for // 2 — 0.45 eV moves up 
to hw = 25.2 meV for /i 2 = 0.9 eV. Interestingly, from this hybridized GMRA-QW plasmon 
peak shift we find that the plasmon energy is not scaled as ~ ^fjMn which is different from the 
observation reported in Ref.[T8l for the ribbon-only case. However, our calculated graphene 
plasmon peak shift (from 14.0 meV to 19.9 me V) for the ribbon-only case (not shown here) 
fully agrees with their observation [18J. The /t 2 -insensitive plasmon peak at hw = 10.7meV 
is associated with the SP excitation at the photon energy hw = a/ e&/ (1 + i&) h£l p \ ~ lOmeV 
and will be addressed in Fig.[6] below. Additionally, when /i 2 = 0.9 eV, we see a new strong 
graphene-like plasmon peak show up at hw = 23.4 meV because one relatively-flattened 
portion [T3"l IT5] on the plasmon dispersion curve at a relative large \q n \ value becomes free 
of Landau-damping. 

To further elucidate the importance of the EM coupling between the graphene plasmon 
in GMRA, as shown in Fig.|2j and the 2DEG plasmon in a QW, we displayed in Fig.|4] 
comparisons for J- p (w) [in (a)] and for the far-field reflectivity spectra TZ p (w), defined by 



Eq. (61 ), [in (b)] as functions of hw by choosing four different array periods: d — 1 /im (blue 
dash-dot-dotted curves), d = 2/zm (black dashed curves), d = 4/xm (red solid curves) and 
d — 8 /im (green dash-dotted curves). Here, we fixed the ratio of the ribbon width W to the 
array period d by C = W /d = 0.5. As indicated by two downward arrows in Fig.[l|a) for the 
shift of one of the two graphene-like plasmon peaks from hw = 19.5 meV for d = 4 /tin to 
hw = 26.6 meV for d = 2 /im, we determine that the hybridized GMRA-QW plasmon energy 
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is approximately proportional to 1/y/d ~ 1/y/W for ( = 0.5. When d = 8 /im, only two 
weakly-split plasmon peaks are visible at fou = 24.1meV and fkv = 24.8 meV, respectively. 
As d decreases from 2 /xm to 1 /xm, the peak at faj = 17.8 meV for d = 2 /xm is pushed 
up to hcu = 24.4 meV. These peak shifts with GMRA period are also reflected in Fig.ffl(b), 
as indicated by the two upward arrows for example, and the steep rise of TZ p (oj) around 
tvjj = lOmeV is associated with the excitation of SP mode in the system. 

In order to clarify the fact that the graphene-like plasmon energy is actually proportional 
to 1/ a/VV, due to the non-locality in the optical response of interacting electrons, instead of 
1/y/d in single-particle modes, we show J~ p {oS) in Fig.[5]for two cases with the same value 
for W but different values for d. For the case with ( = 0.25 and d = 8 /xm, we find a weakly- 
split peak (indicated by two upward arrows), which is aligned with a strong non-split peak 
(indicated by a downward arrow) at £ = 0.5 and d = 4 pan. This clearly demonstrates that 
the graphene-like plasmon peak energy is proportional to 1/ y/\V instead of 1/ y/d, which 
is again in agreement with the observation by Ju et al [18J . The finite width of a ribbon 
introduces a characteristic wavenumber (proportional to 1/W), which enforces a cut-off 
to the short-range coulomb interaction between electrons in a micro-ribbon. A smaller ( 
value implies a weaker coupling between the 2DEG and graphene plasmon excitations. The 
splitting of the plasmon peak for ( = 0.25 is attributed to the contribution of even-integer 
diffraction modes, while only odd-integer diffraction modes can contribute for ( = 0.5, as 



can be verified by the sinc(n'7r£)-terms appearing in Eqs. (44) and (45). 

Figure [6] presents in the left panel the individual and combined effects of SP, 2DEG 
and GMRA as well as the strong EM coupling among these plasmon excitations. The 
periodicity of GMRA opens a minigap between any two adjacent branches (Bloch modes) in 
a folded hybridized-plasmon dispersion curve. These minigaps can be opened either at the 
Brillouin-zone center or at the Brillouin-zone boundary. In the deep sub-wavelength regime, 
only the minigap at the Brillouin-zone center leads to a peak in the far-field transmissivity 
spectrum [ST]. When the doped GaAs substrate is changed into an undoped one (no SP, blue 
dash-dot-dotted curve), only 2DEG and graphene plasmons can exist in this system. In this 
case, there exists strong EM coupling between the graphene and 2DEG plasmon excitations. 
As labeled by the circled numbers, peak-1 and peak-4 are related to the 2DEG-like plasmons, 
while peak-2 and peak-3 are connected to the graphene-like plasmon excitations. Two small 
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kinks on the outer shoulders of peak-2 and peak-3 are a result of two induced anti-crossing 
gaps between coupled 2DEG-like and graphene-like plasmons. If the 2DEG in an InAs 
QW is removed from our system (no QW, black dash-dotted curve), we observe that the 
SP peak is strengthened in comparison with the full system (red solid curve) and shifted 
towards the left (indicated by two solid-line arrows) due to turning-off the 2DEG plasmon 
and its coupling to the SP in the system. Furthermore, the graphene-like plasmon peaks 
superposed on the shoulder of the SP peak become greatly weakened due to the loss of 
the strong coupling between the 2DEG and graphene plasmons. Here, a simple Drude- 
type optical-conductivity model [211 E2] cannot be applied to our system due to the large 
wavenumber involved in the deep sub-wavelength regime for a higher-order diffracted EM 
field. When only the GMRA exists in the system (ribbon only, green dashed curve), a 
downward-shifted dominant plasmon peak at hu = 14.0 meV, along with several kinks (a 
shoulder) below (above) the peak, show up in this panel, as explained in Fig.[2j From 
this figure we know that the graphene-ribbon effect can be best observed from hybridized 
GMRA-QW plasmon modes in the absence of the SP. In the lower panel of Fig.|6j we display 
the calculated transmitted p-polarized E-field intensity |E>(x,^| u)\ 2 at hu = 10.7meV for 
the full system. In this case, a resonant SP effect on the intensity distribution is expected 
around the QW-mediated air/GaAs interface. Indeed, the intensity in the gap between two 
neighboring graphene micro-ribbons, where the QW-mediated air/GaAs interface is hosted, 
decreases greatly due to strong optical absorption by the SP (/r/o = 5meV). At the same 
time, the intensity is built up at two edges of a micro-ribbon and spreads out to the regions 
covered by micro-ribbons. A slight asymmetry in the intensity distribution with respect to 
the left and right edges of the gap region is attributed to a finite incident angle {0i = 30°). 

In the presence of an SP, we compare the transmitted p-polarized E-field intensities 
|E>(x, z\ u)\ 2 in Fig.[7]for the cases of with a QW at hw = 17.4 meV (left panel) or without 
a QW at ftjjj = 16.5 meV (right panel). From the right panel of this figure, we find in 
the near-field regime that a GMRA tends to build up very strong E-field intensities just 
at two edges of a micro-ribbon. On the other hand, the QW would like to spread the E- 
field intensity across the gap region between two neighboring micro-ribbons although two 
close-magnitude maxima in intensity distribution can still be seen between the gap center 
and edges. The competition of these two opposite effects constitutes a strong EM coupling 
between the 2DEG and graphene plasmon excitations, as shown in the upper panel of Fig. [6} 
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In addition, we notice that the E-field intensity for both panels is a little higher at the left 
edge of the gap region for the lower of the two graphene-like plasmon peaks in Fig. [6] 

To find the strongly-diffracted near-field effect in the deep sub-wavelength regime, as 
shown in Fig.[7j on the far-field transmissivity spectrum, we compare in Fig. [8] the calculated 



partial near-field transmissivity spectra F n (q n \ui) for p polarization, defined by Eq. (64), 
with n = ±1 [in (a)] and n = ±2 [in (b)]. It is clear from Figs.[8](a) and (b) that the 
graphene-like plasmon peak at hw = 17.4 me V mainly comes from the contributions of the 
diffraction modes with n = ±1, while the plasmon peak at ftw = 19.5 meV is produced 
jointly by the diffraction modes with both n = ±1 and n = ±2. In general, the peak in 
the near-field transmissivity spectrum (black solid and red dash-dotted curves) can be much 
stronger than that in the far-field spectrum (blue dashed curves) due to the very strong near- 
field intensity as displayed in Fig.[7j The non-locality in the optical response of electrons, 



as given by Eqs. (17) and (19), facilitates the mixing between the EM-field specular and 
diffraction modes, which enables transferring the peak weight from a near field to a far field. 
We also note that the peak strengths for the diffraction modes associated with ±n are close 
in magnitude in this figure. 

Surprisingly, for an inverted structure (QW-ribbon, blue dashed curve), where the 
graphene ribbons become a graphene sheet while the InAs Q W sheet becomes QW rib- 
bons, we find from Fig.[9]that the SP effect is greatly suppressed, compared to the original 
system (G-ribbon, red solid curve), in addition to an overall reduction of peak strength in 
the transmissivity. However, the EM-field reflectivity is found to be enhanced (not shown 
here) in this inverted structure for the range of photon energies shown in the figure. In this 
case, only one very weak peak at huj = 20.1 meV is visible for the inverted structure, which 
is attributed to an order-of-magnitude lower electron density in the Q W compared to that in 
a graphene micro- ribbon. At the same time, the SP- related peak shifts from fvjj = 10.7meV 
to fvjj = 11.7 meV due to switching from the QW-SP coupling to the graphene-sheet-SP 
coupling. 



We display in the upper panel of Fig. 10 the effect of bandgap E G on JF p {ui) with E G = 
(red solid curve), Eq = 0.25 eV (black dash-dotted curve) and Eg = leV (blue dashed 
curve). Since a finite Eq leads to an effective mass (Ea/2v F ) for electrons close to the band 
edge, its effect on ^(uj) is expected to be substantial (see Fig. [9]) as Eq ~ ITwf^f- Indeed, 
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we find from the upper panel that the strengths of the two graphene-like plasmon peaks are 
significantly reduced when E G = leV (for hv F k 2 F = 0.45 eV). Correspondingly, from the 
lower panel of this figure, we see a decrease of near-field intensity as well as an enhanced 
EM coupling between the 2DEG and graphene plasmons, i.e. the E-field intensity spreads 
inward from the two edges to the center of the gap region. 

4. Conclusions 

In conclusion, we have derived in this paper a self-consistent theory involving Maxwell 
equations to determine a total EM field and the density-matrix linear-response theory to 
determine an induced optical-polarization field. We have applied this self-consistent theory 
to study a model system composed of an electronically-isolated two-dimensional electron 
gas and an electron-doped graphene micro-ribbon array on the surface of an n-doped semi- 
infinite semiconductor substrate. The numerical results for the transmissivity and reflectiv- 
ity spectra in the presence of an incident plane wave have been compared and physically 
explained in the deep sub-wavelength regime for various linear-array period, micro-ribbon 
width, doped graphene chemical potential and different configurations of our model sys- 
tem. Our calculations have demonstrated the existence of a tunable EM coupling among 
the localized surface, conventional electron-gas and massless graphene intraband collective 
excitations in the deep sub-wavelength regimes for terahertz frequencies. 
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Appendix A 



Combining Eqs (42)-(49), one can rewrite them into a compact matrix equation as follow: 



M ® u = b 



where the source column vector b is found to be 



(65) 



-A x (q x \ u) 
-A y (q x \ u) 
[(w/c) Xi(Qx, w) + (u/c)(x2(qx, u) 
[(w/c) Xi(g x , w) + (w/c)Cx2(gx, w) 









(ic?7/a;)] A y (q x \u) 
(iw/crj)] A x (q x \u) 



(66) 



If we limit n = 0, ±1, ±2, ■ • • , ±iV, we can write down the unknown column vector u as 



u 



B x (q_ N \ lu) 

B x (q N \u) 
B y (q^ N \u) 

B y (q N \ u) 
C x (q„ N \ u) 

C x (q N \u) 
C y {q^ N \u) 

C y (q N \ u) 



(67) 
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Moreover, elements M(j, f) of the coefficient matrix M introduced in Eq. (65 ) can be written 
out explicitly as 



M(l, j") = I 



1 for f — N + 1 
-1 for f = 5N + 3 
for all other f 



M(2, f) 



1 for f = 3N + 2 
-1 for f = 7N + 4 
for all other j' 



M(3, f) = { 



-% (cp 1)0 /u) - (u/c) \xi(q x , w) + (X2(q x , w)] for j' = 3iV + 2 

-i (cfoflM for J'' = 7iV + 4 

- (u/c) Csinc(m / 7rC) x 2 (qw, w) for j' G [Ji, J 2 ] , [J 2 + 2, J 3 ] 
for all other j' 



where J x = 2N + 2, J 2 = 3iV + 1, J 3 = AN + 2, m' = f - (3N + 2), and 



-i (n 2 a u/cf3 lfi ) - (u/c) \xi(q x , w) + CX2(?*, w)] for / = JV+ 1 



M(4, j") = ^ 



-i (w/c/3 2i o)e s (g ,/32,o, w) 
-{oj/c) Csinc(m'7rC) X2(g m ', 




for j' = 5iV + 3 
for f G [1, N] , [J 1; J 2 ] 
for all other j' 



where Ji = N + 2, J 2 = 2N + 1 and m' = / - (JV + 1). 
For 5 < j < 2N + 4, we have 



M(j, f) = I 



1 for / = j - 4 G [1, iV] or / = j - 3 G [iV + 2, 2N + 1] 

-1 for f = Ji G [4iV + 3, 5iV + 2] or / = J x + 1 G [5iV + 4, 6iV + 3] 
for all other j 1 



with Ji = j + 4iV - 2, while for 2iV + 5 < j < AN + 4, we acquire 
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for f =j - 3 G [2N + 2, 3iV + 1] or f = j - 2 e [3N + 3, AN + 2} 
for f = Ji e [6JV + 4, 7iV + 3] or / = J ± + 1 e [7iV + 5, 8N + 4] 

for all other j' 



with Ji = j + 4iV-l. 

In addition, for AN + 5 < j < 5N + 4, this gives 



-i (cpi,m/u) - (oj/c) [(X2{q m , u) + Xi(Qm, w)] for j' = j - (Ji + 1) 
-i(c/3 2 , m /uj) for j' = j + (Ji - 3) 

- (w/c) Csinc[(m' - m)7rC] X 2 (<7m', w) for j' e [Ji, -h] but j' 7^ J 3 
for all other j 1 



where Ji = 2iV + 2, J 2 = 4iV + 2, J 3 = j - (2JV + 3), m = j-(5iV + 5) and w! = j'-(3N + 2), 
while for 5iV + 5 < j < 6N + 4, it leads to 



-i (cpi, m /oj) - (oj/c) [Cx2{q m , w) + Xi(?m, w)] for j' =j-Ji 

-i(c/3 2 , m /u) for / = j + Ji - 2 

- (w/c) Csinc[(m' - m)7rC] X 2 (<7m', w) for j' G [Ji, J 2 ] but j' 7^ J 3 
for all other f 



where Ji = 2iV + 2, J 2 = 4iV + 2, J 3 = j - (2JV + 2), m = j-(5N + 4) and rri = j'-(3N + 2). 



Finally, for 6N + 5 < j < 7N + 4, this yields 



M(j, f) = I 



-i (n 2 a u/c/3 hm ) - (w/c) l(X2(q m , w) + Xi(?m, w)] for f = j - J2 

-i (Uj/c(3 2 ,m) £s(q m , fam, w) for j' = j - ( Ji + 1) 

- (w/c) Csinc[(m' - m)Tv(] X2(q m >, w) for f G [1, Ji] but j' 7^ J 3 
for all other j' 



where J x = 2N + 1, J 2 = 6N + 4, J 3 = j-(6N + 4), m = j-(7N + 5) and w! = f - (N + 1), 
while for 7N + 5 < j < 8N + 4, one is left with 
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-% {n 2 a u/c(3 hm ) - (u/c) [(X2{q m , w) + Xi(<lm, w)] for / = j - J 2 

. ... -i(u/c(32,m)£s(qm,f32,m,u) for / = j - J x 

J ) = < 

- (w/c) Csmc[(m' - m)ir(] XiiUm'-, w) for j' G [1, -A] but j' 7^ J 3 
for all other f 

where J x = 2N + 1, J 2 = 6N + 3, J 3 =j-(6N + 3), m = j-(7N + A) and w! = j'-(N + l). 



For our numerical calculations in this paper, we take N = 40 to ensure the accuracy of 
the presented results. 
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Fig. 1. (Color online) Schematic representation of a graphene micro-ribbon array (blue) 
with period d and ribbon width W and an InAs quantum well (red). The lower half-space 
(z < 0) is filled with air with refractive index n a — 1, whereas the upper half-space (2 > 0) 
is filled with a doped semi-infinite GaAs bulk having a complex dielectric function e s (q, u) 
or a complex refractive index n s = yje s (q, u). Both the graphene micro-ribbon array [with 
an optical-response function X2{<lx, w)] and the InAs quantum well [with an optical-response 
function Xi(Qx, U} )] sit on the surface [z = 0) of the semi- infinite GaAs bulk. A plane- wave 
electromagnetic field is incident from the z < side with an incident angle 0j. 
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Fig. 2. (Color online) A comparison of the far-field transmissivity spectra T p {oj) in the 
absence of both SP and QW for p polarization with two given linear-array periods: d = 2 fim 
(red solid curve) and d = 4/im (blue dashed curve). Two arrows indicate the shift of two 
corresponding peaks with d. Other parameters in calculations are given in the text. 
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Fig. 3. (Color online) A comparison of the calculated transmissivity spectra Tp^u) for p 
polarization with /i 2 = 0.45 eV (red solid curve) and /i 2 = 0.9 eV (blue dashed curve). Two 
downward solid-line arrows indicate the shift of two corresponding peaks with /z 2 , while one 
upward dashed-line arrow indicates a new peak. Other parameters in calculations are given 
in the text. 




32 



0.5 



0.4 



0.3 



0.2 



0.1 





d=8 um 


1 




1 

s n « ; 

■>\ v /A / \ ' 

x*^ ^ — < V/ \ «' 






■>•"*-->-''*• \ 
'^ N \ \ 

— d=1 um ~ 


i\ ,; \ 

A- / ' 




6=2 um 


(a) 




d-4 um 



10 15 20 25 30 

ft© ( meV ) 



DC 



0.5 



0.4 



0.3 



0.2 



0.1 



it 
\ 

V V 


- d=1 um 

- 6=2 um 

- d=4 uiri^_^_^ 


■ \* 

\\\ 

■ \ . 

\ 

\ / 


— ^ l --T7^a 

x/ r\ v< i 

< " / < ' ■ '. ■' 
\ / \ 1 ' " li ' 


(b) - 


d=8 um t 



10 15 20 25 30 

ft© ( meV ) 

Fig. 4. (Color online) Comparisons of the far- field transmissivity spectra J-" p (o;) [in (a)] and 
the far-field reflectivity spectra TZ p {u:) [in (b)] for p polarization with different linear-array 
periods: d = 1 //m (blue dash-dot-dotted curves), d = 2 //m (black dashed curves), d = 4/xm 
(red solid curves) and d — 8 /im (green dash-dotted curves) . Two arrows indicate the shift 
of two corresponding peaks with d. Other parameters in calculations are given in the text. 
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Fig. 5. (Color online) A comparison of the calculated transmissivity spectra J- p (oS) for p 
polarization with d = 4 /mi and £ = 0.5 (red solid curve) as well as with d = 8 /mi and 
£ = 0.25 (blue dashed curve). The peak indicated by a downward red arrow splits into two 
indicated by two upward blue arrows. Other parameters in calculations are given in the 
text. 
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Fig. 6. (Color online)) A comparison of the calculated transmissivity spectra J- p {uS) is 
presented in the upper panel for p polarization and four different configurations of the 
system, including: (i) with GMRA, QW and SP (full, red solid curve); (ii) Q p i = (no 
SP, blue dash-dot-dotted curve); (iii) Xi{o.xi uj) = (no QW, black dash-dotted curve); (iv) 
Q p i = = Xiilx, w) = (ribbon only, green dashed curve). Two solid-line arrows indicate 
the peak associated with the SP. The circled numbers label four peaks in the figure for the 
case of no SP. For the lower panel, the transmitted p-polarized E-field intensity |E> (x, z\ ui) | 2 
is shown at hu = 10.7 meV for the full system, where the color scale is indicated. Other 
parameters in calculations are given in the text. 
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Fig. 7. (Color online) A comparison of the transmitted p-polarized E-field intensities 
|E>(x, z\ u)\ 2 in the presence of an SP with (left) or without (right) a QW at two indi- 
cated resonant photon energies, where two color scales are given in the left and right panels, 
respectively. Other parameters in calculations are given in the text. 
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Fig. 8. (Color online) Comparisons of the calculated partial near-field transmissivity spectra 
F n (q n \ uj) for p polarization. In (a), we take n = —1 (black solid curve), n = (blue dashed 
curve) and n — 1 (red dash-dotted curve), while in (b) we choose n = —2 (black solid 
curve), n = (blue dashed curve) and n = 2 (red dash-dotted curve). Other parameters in 
calculations are given in the text. 
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Fig. 9. (Color online) A comparison of the calculated p-polarized transmissivity spectra 
Tp{oj) for two compensated structures, including a graphene micro-ribbon array plus an InAs 
quantum-well sheet (G-ribbon, red solid curve) and an InAs quantum-well ribbon array plus 
a graphene sheet (QW-ribbon, blue dashed curve). A blue upward arrow indicates a weak 
peak for the QW-ribbon array. Other parameters in calculations are given in the text. 
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Fig. 10. (Color online) A comparison of the transmissivity spectra T p (oS) is made in the 
upper panel for p polarization with three different bandgaps Eq = (/12 = 450 meV, red 
solid curve), Eq = 0.25 eV (fi2 = 342 meV, black dash-dotted curve) and Eq = leV (fi2 = 
173 meV, blue dashed curve). Two pairs of arrows indicate the shift of a pair of corresponding 
peaks with Eq. In the lower panel, a comparison of the calculated transmitted p-polarized 
E-field intensities |E>(x, z\ u)\ 2 is displayed for Eq = (left) and Eq = leV (right) at 
hw = 19.5 meV and hu = 18.7 meV, respectively. Here, we keep the electron areal density 
(~ k^p) in a graphene micro-ribbon unchanged for different values of Eq, and two color scales 
are indicated in the lower-left and lower-right panels. Other parameters in calculations are 
given in the text. 



39 



